Graphene to Graphane: The Role of H Frustration in Lattice Contraction 



ON 

o 
o 

(N 



S. B. Legoas 1 , P. A. S. Autreto 2 , M. Z. S. Flores 2 , and D. S. Galvao^B 

1 Centro de Ciencias e Tecnologia, Universidade Federal de Roraima, 69304-000, Boa Vista, Roraima, Brazil and 
2 Instituto de Fisica "Gleb Wataghin", Universidade Estadual de Campinas, 
Unicamp, CP. 6165, 13083-970, Campinas, Sao Paulo, Brazil 
(Dated: February, 16th, 2009) 

Graphane is a two-dimensional system consisting of a single planar layer of fully saturated (sp 3 
hybridization) carbon atoms with H atoms attached to them in an alternating pattern (up and 
down with relation to the plane defined by the carbon atoms). Stable graphane structures were the- 
oretically predicted to exist some years ago and just experimentally realized through hydrogenation 
of graphene membranes. In this work we have investigated using ab initio and reactive molecular 
dynamics the role of H frustration (breaking the H atoms up and down alternating pattern) in 
graphane-like structures. Our results show that H frustration significantly contributes to lattice 
contraction. The dynamical aspects of converting graphene to graphane is also addressed. 
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Carbon-based materials are among the most studied 
ones, however the discovery of new structures seems end- 
less, colossal carbon tubes [l[ and graphene Q being re- 
cent examples. 

Graphenes are flat monolayers of carbon atoms in a 
sp 2 hybridization. They have very interesting electronic 
and mechanical properties, which make them one of most 
important subject in materials science today 0,0,13. 

It has been theoretically predicted that a related struc- 
ture, called graphane could exist. Graphane consists 
of a single planar layer structure with fully saturated (sp 3 
hybridization) carbon atoms with H attached to them in 
an alternating pattern (up and down with relation to the 
plane defined by the carbon atoms). The two most sta- 
ble conformations are the so-called chair-like (H atoms 
alternating on both sides of the plane) and boat-like (H 
atoms alternating in pairs) Q (Figure 1). A third mem- 
ber of these two-dimensional planar carbon structures, 
called graphynes @, 0, @|, have also been predicted to 
exist but up to now only molecular fragments have been 
synthesized [7|]. 

Indirect experimental evidences of graphane-like struc- 
tures have been reported 0, [1(3]. More recently, in 



a series of very elegant experiments, Elias et al. 11 1 
demonstrated the existence of graphane formation from 
graphene through hydrogenation with cold plasma. They 
also demonstrated that the process is reversible. These 
fundamental discoveries open new and important per- 
spectives to the use of graphene-based devices since the 
electronic gap values in graphanes could be controlled by 
the degree of hydrogenation 11 , . 

The Elias et al. experiments consisted in exposing 
graphene membranes to H + from cold plasma. The H 
incorporation into the membranes results in altering the 
C sp 2 hybridizations to sp 3 ones. The experiments were 
also done with the membranes over Si02 substrates (only 
one membrane side exposed to H + attacks) and produced 
a material with different properties. 

One very interesting result from their TEM (trans- 
mission electron microscopy) data is the distribution of 



the lattice spacing value of the hydrogenated membranes. 
For an ideal graphane structure theoretical calculations 
indicate a lattice spacing value larger than the graphene 
one @|. Although the TEM data (Fig. 3 from [11]) show 
events with values larger than the graphene one, the ma- 
jority of counts show smaller ones. 

Elias et al. discuss in details a series of experimental 
and material conditions that could explain the contracted 
lattice values, but from the experimental point of view it 




FIG. 1: Structural models, (a) graphene; (b) graphane- 
boatlike; (c) graphane-chairlike. See text for discussions. 
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is very difficult to identify and/or isolate specific contri- 
butions to this effect. Of particular interest it is the role 
of the ordered and/or disordered arrays of H incorpora- 
tion into the membranes. 

In this work we have investigated using ab initio quan- 
tum and classical binding energy bond order (BEBO) 
molecular dynamics methods how the graphene hydro- 
genation is reflected into lattice deformations. 

We have carried out ab initio total energy calculations 
in the framework of the density functional theory (DFT), 
as implemented in the DMol 3 code [131 ] ■ DMol 3 is con- 
sidered state of the art DFT methodology. Exchange 
and correlation terms were treated within the general- 
ized gradient (GGA) functional by Perdew, Burke, and 
Ernzerhof [13 ]. Core electrons were treated in a non- 
relativistic all electron implementation of the potential. 
A double numerical quality basis set with polarization 
function (DNP) were considered, with a real space cut- 
off of 3.7 A . The tolerances of energy, gradient, and 
displacement convergence were 0.00027 eV, 0.054 eV/A 
and 0.005 A respectively. 

We have considered finite and infinite (cyclic boundary 
conditions(CBC)) structures. As we need to investigate 
many structures of different sizes and at different temper- 
atures (molecular dynamics (MD) simulations) the exten- 
sive use of ab initio methods is computationally cost pro- 
hibitive. Thus, for the larger structures and for MD sim- 
ulations, we opted to use a BEBO method like ReaxFF 
H Hill- ReaxFF is a reactive force field developed by 
Adri van Duin, William Goddard III and co-workers for 
use in MD simulations. This empirical method allows the 
simulation of many types of chemical reactions, including 
bond dissociation (making/breaking bonds). 

ReaxFF is similar to standard non-reactive force fields, 
like MM3 [Hj], where the system energy is divided into 
partial energy contributions associated with, amongst 
others, valence angle bending, bond stretching, and non- 
bonded van der Waals and Coulomb interactions 15, lij ]. 
However, one main difference it is that ReaxFF can han- 
dle bond formation and dissociation using bond order 
approach. ReaxFF was parameterized against DFT cal- 
culations. The average deviation between the ReaxFF 
predicted heat of formation values and the experimental 
ones are of 2.8 and 2.9 Kcal/mol for non-conjugated and 
conjugated systems, respectively 

We have carried out geometry optimizations using gra- 
dient conjugated techniques (stopping condition, gradi- 
ent values less than 10" 3 ). For MD simulations we used 
scaling methods, like Beredsen thermostat in order 
to control temperature. 

We started carrying out DMol 3 calculations for the 
infinite (CBC) structures shown in Fig. 1. The cell 
unit parameters were optimized using the Murnaghan 
[20I ] procedure. Sofo, Chaudhari, and Baker [j| in their 
pioneering graphane work considered compacted (inter- 
acting layers) structures. Here, in order to mimic the 



experimental conditions ll| , we have considered isolated 
(non-interacting) layers. The results are displayed in ta- 
ble I. The most stable conformation (chair) has just one 
type of C-C bonds, while the boat conformation has two 
types (carbon bonded to H on opposite side of plane of 
symmetry and those connected to H on the same side of 
it), with their associated H bonds. 



TABLE I: DMol 3 results for the crystalline structures shown 
in Fig. 1. See text for discussions. 

System Energy (Ha) Lattice(A) C-C (A) 



Graphene 


-304.68 


2.465 


1.423 


G-chair 


-309.41 


2.540 


1.537 


G-boat 


-309.38 


2.592 


1.581 






2.509 


1.537 



As we can see from Table I, the lattice parameter values 
for the graphane structures are larger than the graphene 
one. We have also considered the case of the minimum 
unit cell with H atoms parallelly aligned (just one side 
of the membrane). Our results show that the system is 
unstable with the tendency of H 2 recombination and/or 
C-C breaking bonds. These results are in good agreement 
with previous work Q and with the available experimen- 
tal data @,[llj]. 

Considering the stochastic nature of the experiments, 
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FIG. 2: Scheme of the H frustrated domains, 
discussions. 



See text for 
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there is a significant probability of the existence of H 
frustration (Fig. 2). H frustration is a configuration 
where the sequence of alternating up and down H atoms 
is broken (frustrated). This is similar to spin frustration 
in magnetic materials [211 ]. 

In Fig. 2a we show a domain of up and down H atoms. 
From energetic (stability) point of view, after the first 
(up or down) H is incorporated into the C layer, the next 
favorable site is its first inverse neighbor (down or up), 
and so on. If the system is large enough uncorrelated 
domains might be formed (Fig. 2b). As the H coverage 
is continued it could occur that it is no longer possible 
the alternating sequence of up and down H atoms (Fig. 
2c). 

We have investigated finite fragments with and without 
H frustrations. We analyzed the associated geometrical 
changes (Fig. 3) in order to determine whether the lattice 
expands or contracts with relation to the ideal graphane 
structure (table I). We have carried out DMol 3 (table II) 
and ReaxFF (for larger fragments) calculations (table III, 
Fig. 4). 
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FIG. 3: Schematic draw of a graphene fragment before hydro- 
genation. The letters are the reference points for the distances 
displayed in Tables II and III. 

DMol 3 and ReaxFF show similar and consistent re- 
sults. The H frustration produce out of plane distor- 
tions which induce in-plane shrinkage. The net result 
is a decrease of the lattice parameter in relation to the 
ideal graphane value. This effect is amplified when first- 
neighbor H atoms are parallelly aligned (Fig. 4). 

In order to further test these ideas we have also carried 
out ReaxFF molecular dynamics simulations of graphene 
hydrogenation on very large membranes, at different tem- 
peratures and at different H atom densities. 

In Fig. 5 we show representative snapshots from the 
early stages of a simulation at 300 K. The results show 
that significant percentage of uncorrelated H frustrated 
domains are formed in the early stages of the hydrogena- 



TABLE II: Distances between reference points depicted in 
Fig. 3. Frust-1 and Frust-0 refer to parallel (see Fig. 4a) and 
missing (see Fig. 4b) hydrogen atoms in frustrated domains, 
respectively. 



System 


d.A-B 


ds'-c 




du'-A 


dA-C 




Graphene 


9.804 


9.799 


9.804 


9.799 


15.569 


9.974 


G-chair 


9.861 


9.841 


9.882 


9.847 


15.636 


10.050 


G-boat 


9.852 


9.818 


9.852 


9.818 


15.622 


9.977 


frust-0 


9.788 


9.857 


9.876 


9.823 


15.588 


10.003 


frust- 1 


9.740 


9.802 


9.866 


9.786 


15.465 


9.990 



tion process leading to lattice decreased values and exten- 
sive membrane corrugations. These results also suggest 
that large domains of perfect graphane- like structures are 
unlikely to be formed, H frustrated domains are always 
present. The number of these domains seems to be sen- 
sitive to small variations of temperatures and H gas den- 
sities. This can perhaps explain the significant broad 
lattice par ameter distribution values experimentally ob- 
served [ill ]. 

We run annealing cycle simulations to verify whether 
thermal healing occurs. Our preliminary results show 
that, once formed, H frustrated domain are resilient to 
thermal annealing. Further studies are necessary to clar- 
ify this issue. 




FIG. 4: Examples of the optimized geometries for some of the 
frustrated structures listed in tables II and III. (a) frust-1, H 
frustration with parallel first-neighbor H atoms; (b) frust-0, 
H frustration with 'missing' first-neighbor H atoms. Atoms 
in the defect region are shown in ball and stick rendering. 
For view clarity the H atoms outside this region were made 
transparent. 
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TABLE III: Distances between the reference points depicted 
in Fig. 3 calculated with ReaxFF. The number in parenthesis 
indicate the number of frustrated domains in the structure. 
See text for discussions. 



helpful assistance with the ReaxFF code and for his kind 
hospitality in Pennsylvania. 



System 


d.4-B 


d S '_c 


dc-D 


d D i- A 


dA-C 




Graph. 


22.26 


22.26 


22.26 


22.26 


37.08 


22.47 


G-chair 


22.95 


22.94 


22.95 


22.94 


38.40 


23.01 


Frust-1 (13) 


22.21 


22.92 


22.80 


22.90 


37.95 


22.47 


Frust-0 (13) 


22.38 


22.95 


22.77 


22.82 


38.12 


22.37 
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FIG. 5: (a) Representative snapshot of the early hydrogena- 
tion stages from ReaxFF molecular simulations at 300 K. Free 
atomic H atoms are indicated in white and C-bonded ones in 
green, (b) zoomed region indicating H frustrated domains 
formed. The triangle path shows that a sequence of up and 
down H atoms is no longer possible. 



A detailed MD analysis including the effect of many 
layers, substrates and association of inert gas with H 
atoms will be published elsewhere [22l |. 
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